Modulations of the local pairing interaction near magnetic 
impurities in d-wave superconductors 



A. T. R0mer 1 *, S. Graser 2 , P. J. Hirschfeld 3 , and B. M. Andersen 1 
1 Niels Bohr Institute, University of Copenhagen, 
Universitetsparken 5, DK-2100 Copenhagen, Denmark, 
* ar@fys.ku.dk, Tlf: +45 35320451, Fax: +45 35320460 
2 Theoretical Physics III, Center for Electronic Correlations and Magnetism, 
Institute of Physics, University of Augsburg, D-86135 Augsburg, Germany 
^Department of Physics, University of Florida, Gainesville, Florida 32611, USA 

(Dated: March 7, 2013) 

The spin-fluctuation based pairing mechanism has proven successful in explaining 
the pairing symmetries due to Fermi surface nesting of both cuprates and iron-based 
materials. In this work, we study signatures of a spin-fluctuation mediated pairing at 
the local scale. Specifically, we focus on magnetic impurities and calculate both the 
local antiferromagnetism and the resulting modulated pairing interaction. The latter 
gives rise to distinct local enhancements of the superconducting gap in the immediate 
vicinity of the impurities. Our results show that Coulomb- driven pairing naturally 
leads to unusual superconducting gap modulations near disorder potentials. 
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INTRODUCTION 



In the cuprate E^S^CaC^Os+tf a correlation between oxygen dopant positions and local 
enhancements of the superconducting gap has been observed by STM measurements [TJ [2] . 
This correlation was explained at a phenomenological level by introducing a spatially mod- 
ulated pairing strength which is enhanced at the impurity sites J3]. More generally, the 
idea of a spatially modulated pairing interaction caused by the combination of disorder 
and electronic correlations have been somewhat successful in explaining a series of experi- 
ments [3H7]. The microscopic origin of the modulated pairing potential remains, however, 
unsettled despite several recent discussions [8HT3]. We have recently shown that such lo- 
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cal enhancements of the pairing potential and hence the observed superconducting gap is 
a natural consequence of spin-fluctuation mediated pairing [14]. For an overview of spin- 
fluctuation mediated pairing in cuprates and iron-pnictides we refer to Refs. [15j [161 Here, 
we focus on the role of magnetic impurities, and study their effects on the pairing potential 
in rf-wave superconductors. We utilize the recently developed real-space formulation of the 
effective pairing interaction in real-space. Similar to the situation of non-magnetic disorder 
[14], it is found that magnetic impurities may lead to significant local enhancements of the 
superconducting gap. 



MODEL 



In the initial step of the calculation of the effective real-space pairing potential, we obtain 
the spin-resolved charge densities calculated self-consistently in the normal state using a 
mean-field approximation to the one-band Hubbard model 

#o = ^2uAa c JcT + ^2(U(n i(T ) - fin* + J2 V ™p 5 ( t * ~ r *im P )^- (!) 

i,J,(7 1(7 1(7 

Here, c\ a refers to creation of an electron with spin a at lattice site z, and n ia is the number 
operator of spin a particles at site i. Note that this Hamiltonian also contains the impurity 
potential VJ mp at a site placed at position r iimp . A brute force diagonalization of Eq. ([I]) allows 
us to obtain the effective interaction V e ff(i,j) which is due to longitudinal and transverse 
spin fluctuations. Using the approach of Berk and Schrieffer [T7], we obtain a real-space 
formulation of the effective pairing interaction which can be written as [14J 

v eS (ij) = u+ J^mi + . (2 ) 

i - U 2 xnxn to) 1 - u Xit to) 

The susceptibilities entering Eq.Q are real-space matrices given by 

era'/ n\ J \E ma ) — j{Ek, a ') / q n 

Xij {u = 0) = 2^u mia u mja u kj(T ,u kia <— — — — — — — , (3) 

in terms of the eigenvectors u ma and eigenvalues E ma obtained in the diagonalization of 
Eq.Q. 

After the calculation of the effective spin-fluctuation mediated pairing, the density and 
superconducting gap values are calculated self-consistently using a mean-field approach for 
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both density and pairing channels 

where A^- = Ve& ^ (cj^c^ — c^c^). The factor 1/2 arises from the restriction to the singlet 
pairing channel. 

In this real-space study, we investigate the impurity effects on the superconducting gap 
which arise due to a spin-fluctuation mediated pairing mechanism [14] . The self-consistent 
approach allows a detailed study of the co-existence of spin density variations and super- 
conducting gap modulations. We will focus on the effect of point-like magnetic impurities. 



RESULTS AND DISCUSSION 



Prior to the calculation of the effective pairing interaction, the spin densities are cal- 
culated self-consistently for a system containing a single point-like magnetic impurity. In 
Fig. [I] (a-d) the total electron density pi = (n^) + (n^) is shown as a function of lattice 
site. Local density variations occur close to the impurity site, and a small increased charge 
density is evident at the impurity site caused by a finite magnetization. In Fig. [j] (e-h) 
we show the induced magnetization, mi = (n^) — (n^), as a function of lattice site. Local 
antiferromagnetism is induced around the impurity, and is most extended because of smaller 
magnetization induced by weaker impurity potentials. 

The effective superconducting pairing interaction is shown in Fig. [j] (i-1). The interaction 
is calculated in real-space from the expression given in Eq. §2§. Interestingly, the effect of 
the impurity is not only to suppress the interaction at the impurity site. In fact, we see 
an enhancement of the pairing interaction at sites around the impurity. The enhancement 
effect is related to a re-distribution of spin densities which takes place in the presence of 
an impurity as reflected in the total density modulations and local magnetization [14] . The 
effect is most pronounced in the two strongest magnetic impurity cases, see Fig. [j] (k,l). 
In the limit of very strong magnetic impurities, see Fig. [jjh), local antiferromagnetism is 
essentially confined to the impurity site, and still we find significant local gap enhancement 
as seen from Fig. [l](l). It is also known that non-magnetic impurities cause a local enhance- 
ment of the pairing interaction even though no local antiferromagnetism is induced in the 
normal phase of these systems [14] . It therefore turns out that the enhancement effect is not 
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FIG. 1: (a-d) Local charge density of the normal phase. An impurity is positioned at site = 
(13, 13) and the impurity strengths are in units of the nearest neighbor hopping constant t: (a) 
Vimp — 0.5, (b) Vimp = 1 3 (c) Vi mp = 10 and (d) Vi mp = 100. For the results presented here the 
system size is 24 x 24 and the parameters are: U = 2.2, t' — —0.3 in units of t, and the doping 
is x = 0.15. In panels (e-h) we show the magnetization, and in (i-1) the local pairing interaction 
between nearest neighbors, for the same parameters as in (a-d). 

dependent on induced antiferromagnetism and it is local variations of spin densities rather 
than a difference between spin densities that causes the effect. 

Including the effective pairing interaction shown in Fig. [I] (i-1), we calculate the super- 
conducting gap self-consistently by diagonalization of the mean-field Hamiltonian given in 
Eq. Q. As apparent from Fig.[2](a-h), the total density arrangement around the impurity is 
roughly unchanged while more antiferromagnetism is induced around the strongest impurity. 
The superconducting rf-wave order parameter is locally enhanced, see Fig. [2] (i-1), and the 
real-space structure resembles that of the local pairing potential. The enhancement of the 
superconducting gap thus occurs at the borders of regions with spin density modulations 
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FIG. 2: (a-d) Local charge density in the superconducting phase. An impurity is positioned at site 
r km P ~ (13, 13) and the impurity strengths are in units of the nearest neighbor hopping constant 
t: (a) Vi mp = 0.5, (b) Vi mp = 1, (c) Vi mp = 10 and (d): Vimp = 100. For the results presented 
here the system size is 24 x 24 and the parameters are: U = 2.2, t' — —0.3 in units of t and 
doping is x = 0.15. (e-h) Local magnetization for the same impurities and parameters as (a-d). 
(i-1) Magnitude of the local d-wave superconducting order parameter A(i) = ^[A x (i) — A y (i)] for 
the same impurities and parameters as (a-d). 

and the enhancement effect competes locally with antiferromagnetic order. It appears that 
there exists some optimal intermediate impurity strength at which the gap enhancement 
becomes strongest, Fig. [2] (k). 

In Fig. [3] the maximum value of the superconducting gap is shown as a function of 
Coulomb interaction strength, C7, for a magnetic point-like impurity (left) and a non- 
magnetic point-like impurity (right) [14J. In the case of a non-magnetic impurity, no local 
antiferromagnetism is induced in the vicinity of the impurity. However, the enhancement 
effect is still present and the real-space structure of the enhancement resembles that of the 
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FIG. 3: Left: Maximum superconducting gap value as a function of Coulomb interaction strength, 
J7, for point-like magnetic impurities of strengths Vi mp = t, lOt. Also the gap value in the homoge- 
neous system is shown for reference. Right: Maximum superconducting gap value as a function of 
Coulomb interaction strength, J7, for point-like non-magnetic impurities of strengths Vi mp = t, lOt. 

magnetic impurity, shown in Fig. [2] (i-1). One interesting difference occurs for the weak im- 
purity, where a magnetic impurity leads to a stronger enhancement than the non-magnetic 
impurity, compare the blue data points in Fig. |3j 

In Fig. [3] the effect of tuning U becomes apparent. Even though small changes in U will not 
significantly alter the spin density at each lattice site it has a great impact in the RPA-like 
pairing interaction. By increasing the Coulomb interaction we approach the singularities 
in Eq. ^ and therefore the " Stoner enhancement" of the superconducting gap is tuned by 
U. Thus, the Coulomb interaction governs both the local gap enhancements and the de- 
velopment of antiferromagnetism, and the final superconducting gap results from a balance 
between them. 



CONCLUSIONS 



We have shown that a local enhancement of the superconducting gap occurs in the vicinity 
of a magnetic impurity if the pairing interaction is mediated by spin-fluctuations. The real- 
space structure of the gap enhancement is robust to the impurity strength and resembles that 
of a point-like non-magnetic impurity. The effect arises due to local variations in the spin 
densities caused by the impurity. Due to these variations, the spin susceptibilities become 
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inhomogeneous thereby enabling local enhancements of the effective pairing interaction. The 
effect is apparent already from the bare susceptibility|f4|, but is enhanced by correlations 
within RPA. 
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